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We consider the scaling limits for a one-dimensional ran- 
dom growth model, the weakly asymmetric single step Solid- 
on-Solid process. We show that the fluctuation fleld, if con- 
sidered in an appropriate (long) space-time scale, solves the 
Kardar-Parisi-Zhang equation. 
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The problem of determining the macroscopic shape of 
growing processes is a central topic in non-equilibrium 
Statistical Mechanics. Most theoretical progresses have 
been made by analyzing oversimplified models of the 
growth mechanism, which nonetheless catch some of its 
essential physical aspects. A prototypical example of this 
approach is the so-called Eden model A deeper un- 
derstanding of the growth process involves the analysis of 
its deviations from the most probable shape. These fluc- 
tuations phenomena take place on a finer scale, which we 
call here mesoscopic. Accordingly, the randomness in the 
microscopic dynamics is not completely averaged out by 
the coarse-graining procedure, and manifests itself as a 
random term in the mesoscopic equation describing the 
evolution of the fluctuations. 

The long scale behaviour (in mesoscopic units) of the 
fluctuations is believed to be universal, i.e. independent, 
within a suitable class, on the particular growth model. 
The celebrated Kardar-Parisi-Zhang (KPZ) equation |^] 
has been proposed exactly to describe this behaviour. 
A striking feature of growth processes, which has been 
noted already by Eden, is the roughness of the cluster 
surface; the so-called kinetic roughening. At the meso- 
scopic scale this is reflected by the presence, as compared 
to equilibrium phenomena, of large and non-Gaussian 
fluctuations. The KPZ model includes this effect by in- 
troducing a non-linear term in the dynamical equation. 

Having nothing to add on the much debated questions 
of the critical indices and the upper critical dimension 
in the KPZ model, we address here a different, more ba- 
sic, problem. We show how the KPZ equation, which 
is a mesoscopic (continuous) model, can be obtained by 
considering the scaling limit for the fluctuation field of a 
truly microscopic (discrete) random growth model. The 
presence of non-Gaussian fluctuations, as it is clear in 
the theory of critical phenomena, is a highly non trivial 
point. In this contest, the discrete model exhibits in- 
deed Gaussian fluctuations at the hydrodynamical scale, 
whereas they became non-Gaussian on an appropriate, 
longer space-time scale (in microscopic units) . Although 



the discrete model is very peculiar this result confirms 
the universality assumption in the KPZ model. In the 
analysis of the microscopic model we find that the limit- 
ing field solves the KPZ equation with a Wick renormal- 
ization of the non linearity Our discussion will also 
clarify the small gradient hypotheses which is behind any 
derivation of KPZ equation. We stress these results are 
obtained only in the one-dimensional setting, i.e. for in- 
terfaces embedded in R^. The discussion is here kept 
at the physical level of rigour, the more mathematical 
inclined reader is referred to jj] for detailed proofs. 

The microscopic model we analyze is the so-called Sin- 
gle Step Solid-On-Solid process It is an effective 
model for the separation line (interface) between two 
phases. It represents the interface, in a local coordinate 
system, as a single valued function C : Z i— s- Z, where 
Z is the one-dimensional lattice and (^{x) is the inter- 
face height at the point x G Z. The single step require- 
ment is imposed by restricting the allowed configurations 
to those satisfying the constraint \C{x + 1) — C{^)\ = 1- 
We introduce a weakly asymmetric local random dynam- 
ics by the following growth rules: local minima become 
maxima, Ct{^) Ct(^) + 2, with rate 1/2 -f e whereas 
local maxima become minima, Ct(^) Ct(^) ~ 2, with 
rate 1/2. This occurs independently at each site and no 
other transition is allowed. The single step constraint 
is then preserved by the time evolution. This process 
describes a local evaporation/deposition and introduces, 
for e > a growth direction. We also note this process 
can be obtained from the Metropolis dynamics for a two- 
dimensional Ising model in the limit where the inverse 
temperature /3 —> oo and the (positive) external mag- 
netic field H ^ keeping constant the product PH (x e, 
see 1^. In this limit the overhangs are not allowed and 
the interface layer is sharp. In order to obtain the single 
step dynamics also the initial configuration has to be cho- 
sen accordingly, e.g. one can take the plus phase below 
the diagonal and the minus phase above. 

The above dynamics can be represented in terms of the 
Weakly Asymmetric Exclusion Process (WASEP). Let us 
consider a lattice gas with hard core interaction (at most 
one particle per site) and let the particles jump to the 
left (resp. right) with rate 1/2 -|- e (resp. 1/2). We 
denote by 77 = {vi^) , a; G Z} the particles configura- 
tion, 77(2;) G {0, 1}. To WASEP is naturally associated 
a Markov semigroup e*^^ which acts on the functions on 
the configuration space = {0, 1}'^. The kernel of the 
semigroup, e*^= (77, 77'), gives then the transition probabil- 
ity for the transition 77 1-^ 7/. The process is however best 
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characterized by its inGnitesimal generator L^, given by 
L, = Lo + eL- (1) 
where Lq = [L^ + L^)/2 is the symmetric part and 

are the generators of the totaUy asymmetric dynamics. 
In the above equation ry^'^ represents the configuration 
obtained from rj by letting a particle jump from x to y 
and / is a local function on Q. Eq. is simply a formal 
way of writing that in a small time interval dt a particle 
jumps from x to a; ± 1 with probability dt provided there 
is a particle in x and the site a: ± 1 is empty. It is now a 
simple check to verify that the interface dynamics above 
described can be obtained from WASEP by setting [Q 

Ctix) EE J2 i'^vtiy) - 1] (3) 

y<x 

We now introduce a macroscopic coordinate r = ex 
such that the strength of asymmetry coincides with the 
scale parameter. In this coordinate system the interface 
position is given by 

r<(r)=£C.-,(e-V) (4) 

where Cix) is defined by linear interpolation for non inte- 
ger X. Above we scaled the microscopic coordinate x, the 
interface height ( and the microscopic time t diffusively, 
i.e. a; ~ C ~ ^ ~ The results in |^ imply then 
the hydrodynamical limit Assuming that rriQ^r) con- 
verges, as e — > 0, to a continuous function mo(r), then 
TO^(r) mT^r) which solves 

drmr = ^Am, + ^ [1 - (Vm,)2] (5) 

where A (resp. V) is the Laplacian (resp. the gradient). 
Eq. describes the interface displacement as given by 
three factors. A diffusion term which takes into account 
the local relaxation, a constant growth, and a non-linear 
term refiecting the single step constraint: no growth is 
allowed where |Vm(r)| = 1. We also note that, since 
|Vtoo| < 1, Eq. (H) imphes |Vr7ii-| < 1 for all r > 0, as 
it is obvious in the discrete model. 

The hydrodynamical limit is to be interpreted as a Law 
of Large Numbers: the microscopic dynamics is random, 
but in the scaling limit e — s- the deterministic equation 
is obtained. Associated with it there is a Central 
Limit Theorem (CLT) which is stated in the following 
form. Introduce the interface fiuctuations as 

<(r) = [Ce-r(e-V) - E Ce-r(e-V)] (6) 

where E denotes the expectation value and the normal- 
ization y/e is the usual one in CLT. Under some assump- 
tion on the initial configuration (basically the conver- 
gence holds at time 0), converges to the solution of the 
linear Stochastic Partial Differential Equation (SPDE) 



drUr = 2^"^ ~ Vrrir Vut + 1/1 — (Vm-r)^ Wr (7) 

where Wr{r) is the space-time white noise, i.e. 

E {wr (r) Wr' (r')} = S{t - t') S{r - r') (8) 

We stress that in the function to,- is given and it is 
precisely the solution of (|^). The form of equation (Q) 
is easily motivated. The deterministic part is in fact the 
linearization of the hydrodynamical equation (jsj) around 
the trajectory mrir) whereas the prefactor of the noise 
is again due to the single step constraint. 

The results in |^ thus rule out the possibility of ob- 
serving non-Gaussian fluctuations, as predicted by the 
KPZ equation, at the hydrodynamical scale. We argue 
that this will be instead the case when fluctuations on 
longer scale are considered. This possibility has been al- 
ready successfully pursued for some spin-flip (Glauber), 
one-dimensional, long-range (mean-field type) models at 
the critical temperature JlOj . 

In the discussion on the hydrodynamical limit and the 
Gaussian fiuctuations, the condition on the initial config- 
uration were rather general; for instance one can choose 
a random 770 distributed with a slowly-varying product 
measure such that 

E{rjo{x)} = poiex) (9) 

and then take Co(a^) as given by (|^). In this case one ob- 
tains mo such that Vrng — 2po~l and a Gaussian process 
for uq; its variance is determined by pq. In order to derive 
the KPZ equation we have instead to prepare more care- 
fully the initial state We assume Co{x) to be a small 
perturbation of a straight line, furthermore the perturba- 
tion is so slowly-varying that it changes only on the scale 
e^^, which is much longer than the hydrodynamical scale 
previously discussed. Denoting by Co(^) = ^ tan 9 the 
straight interface we thus assume there exists an Holder 
continuous function /iq : R i-^ R (it may be random) 
such that 

£ [Co(a^)-Co(a^)] ~ hoie^x) (10) 
where the approximate identity means that for any r G R 
lim e{Co(£-M-Co(£"'0} -/io(r) (11) 

This is the precise formulation, in our contest, of the 
small gradient condition in the KPZ model. 

A naive convergence result will then be the following. 
Assuming (^0|), the interface position at time t is 

(t{x + etta.n0) w Ct(a^ + £tta,n0) + e^^h^nie'^x) (12) 

where Ct the evolution of the straight interface, i.e. 
Ctix) = Coi^) + £t{l — tan^ ^)/2, and hr = hr{r) is the 
solution of the KPZ equation 
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drhr = ]^^hr ~ Kf + V^l - t&T? 9 Wr (13) 

We stress the fluctuation normalization in (|l^), even if 
the non-linear Eq. (^3|) is obtained, is still the CLT [non 
anomalous) normalization. We also note the space-time 
scaling is still diffusive, x ^ e~^,t ~ e~^. As usual ||] 
a Galilean transformation in (^|) is performed to obtain 
the KPZ equation. 

The argument leading to the KPZ equation is the fol- 
lowing. By Eq.s (0), (||) and (||), we have 

LeCix) = i(l +e)V+V-C(a;) + | [l ^ V+C(a;)V-C(x)] 

(14) 

where V^({x) = ±[({x ± 1) — C{^)] are the discrete gra- 
dients and V+V" the discrete Laplacian. Inserting the 
ansatz ( [T^ ) in (|lj) and dropping higher order in e, one 
obtains the deterministic part on the right hand side of 
(p^). Analogously, the noise term can be identified by 
the computation 

= Six-y){l~\/+Cix)V-ax) 

+e [I + V+V'Cix) - V+Cix)V'Cix] } 

The main point in the above formal argument is of 
course the vahdity of ( [l^ ) under the assumption (10), 
which is a propagation of chaos result. For WASEP this 
has been proven in a strong form |^ , but those results do 
not hold in the time scale we are interested in. In fact we 
implicitly assumed h-j-^r) to be a smooth (differentiable) 
function, but a simple analysis of ( [T^ ) fl^ , shows this 
cannot be the case and that r l—^ hr{r) can be at most 
a-Holderian, a < 1/2. 

As emphasized in fl^ , a correct way to define the pro- 
cess hr is through the Cole-Hopf transformation. For 
notation simplicity we discuss next only the case of a flat 
interface, 9 — 0, so that no Galilean transformation in 
( p^ ) is needed. Let V-'r be the solution of the stochastic 
heat equation 



dr^r = 7: AV't + V't Wr 



(15) 



which is linear and (in one space dimension) completely 
meaningful in the space of continuous functions when 
the stochastic differential is interpreted in the Ito sense. 
Furthermore, if the initial condition is positive, then the 
process is strictly positive for any t > as proven in 
[0. We can then define (rigorously) hr ~ Xog^r, which 
(formally) solves the KPZ equation with a Wick renor- 
malization of the non linearity, i.e. 



drhr = ^^/Jr - ^ : (V/l^)^ : +Wt 



(16) 



where the Wick product can be informally written as 
: {Vhrf ■■ = {"s/hrf - ^(0). The extra term (5(0) arises 
because Ito's rules of stochastic calculus are to be used. 

It is a fortunate coincidence that an analogous of the 
Cole-Hopf transformation, first introduced by Gartner 
Q , exists also for WASEP. Given C,t let us define a trans- 
formed process by 



^t{x) = exp{-7eCt(a;) -I- A^i} 



(17) 



where 7^ — log ^/\ + 2e and = 1 + e — ^/l + 2e. By a 
straightforward computation, it is verified that solves 
the stochastic equation 



dUx) = — V+V-6(a;)di + dMt{x) 



(18) 



where the first (linear) term is just the discrete Laplacian 
and the noise term Mt is characterized by 



dt 



{M{x),M{y))t^2e^5{x-y)Uxf 



(19) 



X [t1s{x){1 - ?7,(x + 1)) + e-'T=?7,(a; + 1)(1 - 77,(x))] 

The problem of showing the convergence of the inter- 
face fluctuations to the KPZ equation is thus reduced in 
proving the convergence of S.e-iri^^'^^) to iprir). Since 
the first, linear, term is precisely the discretization of the 
Laplacian in (|l^), it comes for free. To identify also the 
noise, one has to show, using the basic hypotheses (p^), 
that the term in the square bracket on the right hand 
side of (19) converges to 1/2. As rj{x) is a small pertur- 
bation of the product measure with density 1/2 this is 
intuitively clear and it can be rigorously justified 

By taking the inverse Cole-hopf transform we get 



Ct(x)~V,t + e'^h,.t{e'x) 



(20) 



where Vg = \elle = e/2 — £''/4!-t-0(e'*) and hr{r) is now 
associated with the renormalized KPZ equation (p^). 

A few comments are due. We first address the Wick 
renormalization. We note that the product measure with 
a constant density 1/2 is stationary for WASEP, by (||) 

has then a steady state. The expectation of the non 
linear part in ( p^ ) in this steady state is 



Ei/2{V+Ct(x)V-Ct(^)} 
= Ei/2{(27?t(x)-l)(27yt(2; 



1)-1)} = 



which is the same property that defines the Wick prod- 
uct. This shows very clearly that the non linearity in the 
KPZ equation cannot really be simply (Vft-r)^, because 
the latter would have a positive expectation. 

In a more general framework the transition from a 
microscopic description to the macroscopic equation is 
carried out by using some local equilibrium hypotheses 
(proven to hold in number of different models) which al- 
lows to replace microscopic averages (in this case e(^t{x)) 
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by macroscopic quantities {hr{r)). We have here dis- 
cussed a situation where a naive local equilibrium fails, 
but it is sufficient a minimal modification of the macro- 
scopic term, i.e. the Wick product, to get the correct 
answer. In our approach this fact is somehow hidden in 
the Cole-Hopf transformation, a direct argument would 
be much neater and helpful in suggesting how to deal 
with the general cases where no Cole-Hopf transforma- 
tion is available. 

Comparing ( |2C| ) with (|l2|) we see the extra term 
(which is finite in the scaling limit) in the drift veloc- 
ity Ve- To explain it we have to recall the well-known 
connection between the KPZ equation and the directed 
polymers in a random environment. The solution of 
the stochastic heat equation (|l5| ) can in fact be inter- 
preted as the partition function of directed polymers in 
the quenched random environment described by the po- 
tential w. Accordingly, the quenched free energy per unit 
of length is given by |l5| 

hm . 1 (21) 

T-foo r 4! 

We now recall that the solution of KPZ equation has 
been defined by hr — log ijjr ■ Since the expectation of C,t 
in the steady state is given by Ei/2{Ct(2;)} = we 
see that the second term in Ve is indeed needed to get 
the correct result. We can also view the above discussion 
as a non-replica argument (but still not a rigorous proof) 
to get the exact value of the quenched free energy. 

We finally comment on the most unsatisfactory point 
in our approach, that is the fact e is hoth the scale pa- 
rameter and the asymmetry strength. It would be tempt- 
ing to try the asymmetric exclusion process, where the 
asymmetry is kept fixed and independent on e. However 
in this case the correction to the hydrodynamic (which is 
here the inviscid Burgers equation in the Euler scaling) 
are expected to appear at time e~^/^, before the random 
noise in the dynamics shows up Q. In two dimension 
logarithmic corrections are instead expected. In three or 
more space dimension there is a rigorous analysis beyond 
the hydrodynamical scale and the limiting equation 
is deterministic. The interpretation in term of the inter- 
face dynamics is anyhow purely one-dimensional. 

To summarize, we have analyzed a one-dimensional 
growth process beyond its first nontrivial hydrodynam- 
ical scale (which is here the diffusive scale) and found 
a scaling limit where the fluctuations evolve exactly ac- 
cording to the KPZ equation. In the Renormalization 
Group terminology this can be stated saying that the ir- 
relevant term in the microscopic dynamics do not survive 
the limit, thus proving the universality assumption. In a 
more precise formulation, which takes full advantage of 
the Cole-Hopf transformation, the Wick renormalization 
of the non-linearity appears naturally. 
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